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Abstract. We explore the adiabatic pulsational properties of massive white dwarf stars with hydrogen-rich en- 
velopes and oxygen/neon and carbon/oxygen cores. To this end, we compute the cooling of massive white dwarf 
models for both core compositions taking into account the evolutionary history of the progenitor stars and the 
chemical evolution caused by time-dependent element diffusion. In particular, for the oxygen/neon models, we 
adopt the chemical profile resulting from repeated carbon-burning shell flashes expected in very massive white 
dwarf progenitors. For carbon/oxygen white dwarfs we consider the chemical profiles resulting from phase sepa- 
ration upon crystallization. For both compositions we also take into account the effects of crystallization on the 
oscillation eigenmodes. We find that the pulsational properties of oxygen/neon white dwarfs are notably different 
from those made of carbon/oxygen, thus making asteroseismological techniques a promising way to distinguish 
between both types of stars and, hence, to obtain valuable information about their progenitors. 
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1. Introduction 

White dwarfs are the most common end-product of stellar 
_ evolution. Most of the presently observed white dwarfs are 
1 post- Asymptotic Giant Branch (AGB) remnants, the core 
of which consists of the main ashes of helium burning, ba- 
sically a mixture of carbon and oxygen. However, there is 
growing theoretical evidence suggesting that white dwarfs 
with masses larger than ~ 1.05 Mq could have cores com- 
posed primarily of oxygen and neon (Garcia-Berro & Iben 
1994; Ritossa et al. 1996; D'Antona & MazziteUi 1996; 
Garcfa-Berro et al. 1997; Iben et al. 1997). In particular, 
Garcfa-Berro et al. (1997) found that when the core mass 
of the 9 Mq white dwarf progenitor exceeds ^ 1.05 M©, 
and before reaching the thermally pulsing phase at the 
AGB tip, carbon is ignited off-center in semidegener- 
ate conditions. As a result of repeated carbon-burning 
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shell flashes that eventually give rise to carbon exhaus- 
tion in the degenerate core, it is found that at the end 
of the carbon burning phase the star would be left with 
an oxygen/neon core with trace amounts of carbon and 
other heavier chemical species. After considerable mass- 
loss episodes, the progenitor remnant is expected to evolve 
into the central star of a planetary nebula and ultimately 
into a white dwarf with an oxygen/neon core of mass 
~ 1.05 Mq. a possible observational counterpart of these 
massive white dwarfs would be the single massive white 
dwarf LHS 4033, which has a mass of ^ 1.32 Af© (Dahn 
et al. 2004). Other possible massive white dwarfs hosting 
oxygen/neon cores would be the magnetic white dwarf 
PG 1658-^441 (Schmidt et al. 1992; Dupuis et al. 2003) 
— with a mass of ~ 1.31 Mq — and the highly magnetic 
white dwarf RE J0317~853 (Barstow et al. 1995; Ferrario 
et al. 1997), which has a mass of ~ 1.35 Af©. It should 
be noted, however, that there are alternative evolution- 
ary channels that could eventually lead to massive white 
dwarfs with carbon/oxygen cores. These scenarios involve 
the merging of two otherwise light-weight ordinary white 
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dwarfs. It has been shown (Segretain et al. 1997; Guerrero 
et al. 2004) that in this case nuclear reactions are not able 
to modify considerably the composition of the core of the 
remnant. 

Many white dwarfs exhibit multiperiodic luminosity 

variations caused by gravity-modes of low harmonic degree 
{i < 2) and periods ranging from roughly 100 s to approxi- 
mately 1200 s. Over the years, the study of the pulsational 
patterns of variable white dwarfs through astcroscismolog- 
ical techniques has become a very powerful tool for prob- 
ing the inner regions that would be otherwise unacces- 
siblc to direct observations. In particular, for hydrogen- 
rich variable DA white dwarfs — also known as ZZ Ceti 
stars — such technique has proved to be very successful in 
providing independent valuable constraints to their funda- 
mental properties. Examples of these are their core compo- 
sition, the outer layer chemical stratification or the stellar 
mass — sec, for instance, Pfeiffer et al. (1996), and Bradley 
(1998, 2001), amongst others. Very recently, increasing at- 
tention has been paid to study the asteroseismological 
properties of massive ZZ Ceti stars, since it opens the 
interesting possibility of probing the physical mechanisms 
operating in their very dense interiors and, particularly, to 
obtain useful information about the crystallization process 
which occurs in their cores (Montgomery & Winget 1999). 
This has been motivated by the discovery of pulsations in 
the star BPM 37093 (Kanaan et al. 1992), a massive ZZ 
Ceti star which has a stellar mass of ~ 1.05 and an 
effective temperature T^g ~ 11800 K, and which, there- 
fore, should have a sizeable crystallized core (Winget et al. 
1997). It is interesting to note at this point that this mass 
is very close to the theoretical lower limit (^ 1.05 Mq) 
for an oxygen/neon white dwarf to be formed (Salaris et 
al. 1997; Gil-Pons et al. 2003). Indeed, Gil- Pons et al. 
(2003) find that for carbon burning to take place in the 
corresponding progenitor, its mass in the main sequence 
should be ~ 8.1 Mq, leaving a massive oxygen/neon white 
dwarf of ~ 1.05 Mq. Hence, BPM 37093, could be either a 
carbon/oxygen or an oxygen/neon white dwarf, depending 
on the precise value of its mass. 

We do not intend to perform a detailed modelling of 
the pulsational characteristics of BPM 37093, since we 
consider it to be beyond the scope of the paper. Instead, 
the aim of the present paper is to assess the adiabatic 
pulsational properties of massive white dwarfs with car- 
bon/oxygen (CO) and oxygen/neon (ONe) cores. More 
specifically, our main goal is to explore the possibility 
of using white dwarf asteroseismology to distinguish be- 
tween both types of stars. To this end, we compute the 
cooling of massive white dwarf models for both core com- 
positions taking into account the evolutionary history of 
their progenitors and the chemical evolution caused by 
time-dependent element diffusion during the evolution. 
The paper is organized as follows. In Sect. 2, we describe 
the main characteristics of our CO and ONe model white 
dwarfs. Special emphasis is given to the formation of the 
core chemical profile during the pre- white dwarf evolution- 
ary stages. The pulsational predictions for the both types 



of white dwarf sequences are described in §3. Finally, in 
§4 we briefiy summarize our main findings and draw our 
conclusions. 



2. Input physics and evolutionary models 

In this work we compute the evolution and the pulsational 
properties of massive white dwarfs with ONe cores. In pur- 
suing this goal we adopt the chemical profiles obtained by 
Garcia-Berro et al. (1997) — see §2.2 for a detailed dis- 
cussion. Additionally, since a major aim of our work is 
to compare the theoretical pulsational spectra of our ONe 
white dwarf models with those of their CO analogs, we also 
compute the evolution of a white dwarf with a CO core of 
similar mass. The most relevant results of this calculation 
are presented in §2.1 below. An important aspect of the 
present study common to both sequences concerns the evo- 
lution of the chemical abundance distribution caused by 
diffusion once the white dwarfs arc formed. This is partic- 
ularly relevant for the study of the pulsational properties 
of white dwarfs. Indeed, element diffusion turns out to be a 
key ingredient as far as mode trapping in pulsating strati- 
fied white dwarfs is concerned (Corsico et al. 2001; Corsico 
et al. 2002). We have used a time-dependent treatment for 
multicomponent gases (Burgers 1969) which accounts for 
gravitational settling and chemical and thermal diffusion 
for ^H, ^He, ^He, ^^C, and ^^O species, the dominant con- 
stituents outside the core of our white dwarfs models — 
see Althaus et al. (2001a, 2001b) for additional details. In 
this way, the trace element approximation usually invoked 
in most ZZ Ceti studies is avoided. 

2.1. Carbon/oxygen white dwarf models 

Our massive white dwarf models with a CO core have 
been computed with the evolutionary code described at 
length in Althaus et al. (2003). The code is based on an 
up-to-date and detailed physical description. Briefly, the 
code uses OPAL radiative opacities for arbitrary metal- 
licity from Iglesias & Rogers (1996) and from Alexander 
& Ferguson (1994) for the low-temperature regime. The 
equation of state for the low-density regime comprises 
partial ionization for hydrogen and helium composi- 
tions, radiation pressure and ionic contributions. For the 
high-density regime, partially degenerate electrons and 
Coulomb interactions are also considered. Under degen- 
erate conditions we use an updated version of the equa- 
tion of state of Magni & Mazzitelli (1979). Neutrino emis- 
sion rates and high-density conductive opacities are taken 
from the works of Itoh and collaborators — see Althaus 
et al. (2002) . A nuclear network of 34 thermonuclear reac- 
tion rates and 16 isotopes has been considered to describe 
hydrogen (proton-proton chains and CNO bi-cycle) and 
helium burning. Nuclear reaction rates were taken from 
Caughlan & Fowler (1988) except for the ^^C(a, 7)^60 re- 
action rate, for which we adopted that of Angulo et al. 
(1999). 
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Fig. 1. The internal chemical profile corresponding to a 
1.06 CO model in which chemical rehomogenization 
induced by crystallization has been neglected. The gray 
area marks the domain of crystallization. Mc means the 
crystallized mass of the model. 

The treatment of the change of the abundances dur- 
ing the pre-white dwarf evolution is an important as- 
pect to be considered. In particular, our code uses a 
time-dependent scheme for the simultaneous treatment 
of chemical changes caused by nuclear burning and con- 
vective, salt finger and overshoot mixing, which are de- 
scribed as diffusion processes (Althaus et al. 2003). We 
have used the Schwarzschild criterion for convective stabil- 
ity. However, the occurrence of additional mixing beyond 
what is predicted by such criterion is suggested by both 
theoretical and observational evidence. In particular, ex- 
tra mixing episodes (particularly mechanical overshooting 
and/or semiconvection) taking place beyond the formally 
convective boundary towards the end of central helium 
burning have a large infiuence on the carbon and oxygen 
distribution in the core of white dwarfs — see Straniero et 
al. (2003) for a recent discussion. The occurrence of such 
mixing episodes leaves strong signatures on the theoreti- 
cal period spectrum of massive ZZ Ceti stars (Althaus et 
al. 2003). Therefore, we have allowed for some mechanical 
overshooting by following the formalism of Herwig (2000). 
Our mixing scheme allows for a self-consistent treatment 
of diffusive overshooting in the presence of nuclear burn- 
ing. In particular, we have considered exponentially decay- 
ing diffusive overshooting above and below any convective 
region, including the convective core (main sequence and 
central helium burning phases), the external convective 
envelope and the short-lived helium-fiash convection zone 
which develops during the thermal pulses. Finally, con- 
vection is treated according to the mixing length theory 
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Fig. 2. Same as Fig. ^ but in the case in which chemi- 
cal rehomogenization induced by crystallization has been 
accounted for. 



for fiuids with composition gradients (Grossman & Taam 
1996) that applies in the convective, semiconvective, and 
salt finger instability regimes. The mass of the resulting 
CO white dwarf is « 0.94 M©. However, in order to per- 
form a direct comparison with the ONe white dwarf model 
previously mentioned we have artificially scaled the stellar 
mass of the CO model to 1.06 M© — approximately the 
mass of BPM 37093 — but taking into account the core 
chemical distribution expected for a progenitor star with 
a stellar mass correspondingly larger. 

Crystallization in the core of our ZZ Ceti models and 
its effects on their pulsational pattern constitute an cen- 
tral point of this paper. Crystallization sets in when F = 
/rk^T = 180, where r is the radius of the Wigner- 
Seitz sphere. In addition, the chemical redistribution in 
the fluid above the solid core induced by phase separation 
upon crystallization (Garcia-Berro et al. 1988; Segretain 
et al. 1994) has been taken into account by adopting the 
phase diagram of Segretain & Chabrier (1993) and follow- 
ing the procedure described in Salaris et al. (1997) — see, 
also, Montgomery et al. (1999). Nevertheless, for the sake 
of completeness we have also computed a cooling sequence 
in which phase separation has been neglected. 

In our evolutive sequence, the CO white dwarf model 
begins to crystallize at Tcff ~ 16600 K. In Fig. ^we show 
the chemical profile of a CO white dwarf model at the 
ZZ Ceti stage, in which the chemical redistribution due to 
phase separation has been neglected. The effective tem- 
perature of the model is Tcff « 11800 K and their crystal- 
lized mass fraction amounts to 0.63 (the crystallized core 
is shown as a gray zone). Note the presence of a pron- 
nounced step at log(l — Mr/M*) « —0.7 in the oxygen 
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Fig. 3. The internal chemical profile corresponding to our Fig. 4. Same as Fig. El but for an ONe white dwarf model 
initial ONe white dwarf model in terms of the outer mass in the ZZ Ceti instability strip. The crystallized region is 



fraction. 



displayed as a gray area. 



profile. This feature, which reflects the occurrence of over- 
shoot episodes prior to the formation of the white dwarf, 
leaves strong imprints on the theoretical period spectrum 
(Althaus et al. 2003). Note as well that the external chem- 
ical interfaces, including the He-core transition are very 
smooth as consequence of the element diffusion processes 
acting during the evolution. 

Fig. 121 displays the chemical profiles of a CO ZZ Ceti 
model in which chemical redistribution due to phase sepa- 
ration has been considered at Tcs « 11800 K. In this case 
the crystallized mass fraction amounts to 0.55, a slightly 
smaller value than that of the case in which phase sep- 
aration was neglected (Fig. QJ. This is due to the larger 
abundance of carbon in the overlying fluid layers resulting 
from the oxygen enhancement in the crystallized core of 
the model. Note that, as a consequence of the mixing re- 
sulting from phase separation, the overshoot-induced step 
in the innermost oxygen profile has been completely wiped 
out. Consequently, the theoretical period spectrum for this 
model should be much less featured when compared with 
that of the case in which chemical rehomogenization has 
been neglected. This feature will be discussed in depth in 
§3. 

2.2. Oxygen/neon white dwarf models 

The evolutionary stages — from the zero-age main se- 
quence, through helium burning and off-center carbon ig- 
nition and burning in partially degenerate conditions up to 
the thermally pulsing AGB — leading to the formation of 
ONe white dwarfs, as well as the input physics employed 
in the calculations are described in Garcia-Berro et al. 



(1997) and references therein. Here, we restrict ourselves 
to summarize the main characteristics of the models and, 
particularly, their internal composition. When the carbon 
burning phase is already finnished, the star is left with 
a ~ 1.06 core primarily composed of ^^O and ^'^Ne, 
with trace amounts of ^^C, ^'^Na and ^^Mg. Surrounding 
the core there is a buffer composed mainly of a mixture of 
^^C and ^^O and, on top of it, the overlying hydrogen- and 
helium-rich envelope layer. The initial envelope chemical 
profile — which, after diffusion, results into a pure hydro- 
gen layer — is shown in Fig.O We have obtained an initial 
configuration for our ONe white dwarf model by simply 
scaling the internal chemical profiles built-up during the 
pre-whitc dwarf evolution to the structure of the massive, 
hot CO white dwarf model considered previously. In this 
way we obtain a good starting configuration with a mass 
close to that of BPM 37093 and the appropriate chemical 
profile. Because we are interested in the comparison of the 
pulsational properties of highly- evolved CO and ONe ZZ 
Ceti models with the same stellar mass, this procedure is 
enough for our purposes. 

In our evolutive sequence, the ONe white dwarf model 
begins to crystallize at Tcff « 20700 K, long before reach- 
ing the ZZ Ceti instability strip. Fig. 01 shows the result- 
ing chemical structure when the model has reached the 
ZZ Ceti instability strip. At this point of the evolution, 
the percentage of crystallized mass fraction amounts to 
~ 90%. Two important points deserve additional com- 
ments. First, as it is the case for CO white dwarfs, the 
effects of chemical diffusion in the external chemical abun- 
dance distribution is clearly noteworthy. Indeed, the outer 
chemical interfaces are markedly smoothed out by dif- 
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fusion. Second, a rehomogenization process induced by 
Rayleigh- Taylor instabilities has led to a plateau in the 
innermost chemical profile by the time the ZZ Ceti stage 
is reached. Note that as a result, a strong step in the 
chemical profile at log(l — M-^/M^,) k, —1.25 has arisen 
after rehomogenization, a feature expected to be the most 
important ingredient in determining the structure of the 
period pattern of the model. It is as well important to 
mention here that ONe white dwarf models do not expe- 
rience a significant chemical redistribution due to phase 
separation because the charge ratio between Ne and O is 
much smaller than that of O and C and, hence, the phase 
diagrams of Segretain & Chabrier (1993) do not predict a 
sizeable enrichment in Ne in the solid phase. Moreover, the 
evolutionary calculations of Ritossa et al. (1996), Garcia- 
Berro et al. (1997) and Iben et al. (1997) predict that 
after the carbon burning phase and when the thermally 
pulsing phase ensues there is a non-negligible fraction of 
unburnt carbon in the oxygen/neon core. The mass frac- 
tion of unburnt carbon in the ONe core can be as high 
as Xc ~ 0.01 — see FigOl The same happens for the 
case in which a merger of two light-weight white dwarfs is 
involved (Guerrero et al. 2004). Hence, the calculation of 
the effects of Ne sedimentation upon crystallization would 
require a ternary phase diagram (Segretain 1996) which, 
at present, is not yet well known and, moreover, is beyond 
the scope of this paper. 

3. Pulsational calculations 

We have carried out an adiabatic pulsational analysis of 
our white dwarf models with the help of the same Newton- 
Raphson pulsational code employed in Corsico et al. (2001; 
2002) and described in detail in Corsico (2003). This code 
is coupled to the evolutionary code previously described. 
However, for the purposes of this work, we have done the 
appropriate modifications to the pulsational code in or- 
der to properly handle the effects of crystallization on 
the oscillation eigenmodes. Briefly, the boundary condi- 
tions at the stellar center (when crystallization has not 
yet set in) and surface are those given by Osaki & Hansen 
(1973). However, when the core of the white dwarf under- 
goes crystallization we switch the fluid internal boundary 
conditions to the so-called "hard sphere" boundary con- 
ditions (Montgomery & Winget 1999). Within this ap- 
proximation the nonradial eigenfunctions are inhibited to 
propagate in the crystallized region of the core. Following 
previous studies of white dwarf pulsations, the normaliza- 
tion condition adopted is S^-c/r = 1 at the stellar surface, 
being ^i- the radial component of the displacement. The 
oscillation kinetic energy (i?kin) is computed according to 
Eq. 1 of Corsico et al. (2002), whereas the weight function 
is computed as in Kawaler et al. (1985). Finally, we derive 
the asymptotic spacing of periods (AP^) as in Tassoul et 
al. (1990): 
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Fig. 5. The squared Brunt- Vaisala frequency (iV^) — 
solid line — in terms of the outer mass fraction, corre- 
sponding to the same CO white dwarf model shown in 
Fig. n Inset: the Ledoux term B. For the sake of com- 
pleteness, we also show using a dotted line the acoustic 
(Lamb) frequency corresponding to ^ = 2. The gray area 
marks the domain of the crystallized core. Here, because 
of the hard-sphere approximation, the detailed shape of 
the Brunt- Vaisala frequency is irrelevant to pulsations and 
hence it is not shown. 



27r^ 



1 



dr 



(1) 



where N is the Brunt- Vaisala frequency, and ri and r2 
are, respectively, the radii of the inner and outer bound- 
aries of the propagation region of the modes. Hence, when 
the core is already partially crystallized ri coincides with 
the radius of the crystallization front, ri = r{Mc). It is 
important to realize that the turning point ri is no longer 
a fixed value but, instead, it is a function of the crys- 
tallized mass, Mc and, hence, of the temperature of the 
(nearly) isothermal core. As a result, as the white dwarf 
cools down, the internal boundary at ri moves outward, so 
the integral J^^ Ndr/r decreases, and, consequently, the 
period spacing increases and the periods themselves do so 
— see Figs. 7 and 8 of Montgomery & Winget (1999). 



3.1. Brunt-Vaisala frequency 

The Brunt-Vaisala frequency is computed according to the 
procedure of Brassard et al. (1991). This numerical treat- 
ment accounts explicitly for the contribution to N from 
any change in composition in the white dwarf models by 
means of the Ledoux term B. This is an important aspect 
in connection with the phenomenon of mode trapping and 
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Fig. 6. Same as in Fig. [3 but for the CO white dwarf 
model analyzed in Fig. |21 



mode confining (Brassard et al. 1992; Corsico et al. 2002). 
Specifically, the Brunt- Vaisala frequency is given by: 



P Xp 



(2) 



The Ledoux term B, for the case of a multicomponent 
plasma, is given by 



B 



1 '4^' d\nX, 



J2 ^ 



(3) 



Here, p is the pressure, Xt iXp) denotes the partial loga- 
rithmic derivative of the pressure with respect to T (p), 
V and Vad are the actual and adiabatic temperature gra- 
dients, respectively, Xi is the abundance by mass of the 
chemical species i, and 



dlnp 
dlnXi 



(4) 



We begin by examining Fig. [3 where the square of the 
Brunt- Vaisala frequency corresponding to the same CO 
white dwarf model displayed in Fig.^is shown. It is worth 
recalling that in this model phase separation upon crystal- 
lization was disregarded. Note as well that in our models, 
the shape of the outer chemical interfaces is assessed using 
a time-dependent element diffusion scheme. Consequently, 
and not surprisingly, our calculations predict the presence 
of very smooth bumps in the profile of N'^ in the external 
layers of the model (Corsico et al. 2002). Hence, and con- 
trary to what it is usually found when the so-called "trace 
element approximation" is used, in our calculations these 



features play a minor role in the mode trapping and confin- 
ing characteristics of the theoretical pulsational spectrum. 
Instead, it is the chemical structure of the core what de- 
termines the structure of the period spectrum, as already 
shown in Althaus et al. (2003). Specifically, the pulsation 
characteristics of this model are fixed by the presence of 
the pronounced peak in N at log(l — M^/M^,) « —0.7, a 
feature directly related to the overshoot-induced step in 
the chemical profile (see Fig. 

When phase separation upon crystallization is taken 
into account for the CO white dwarf models, chemical re- 
distribution in the overlying fluid layers is expected to 
take place as a consequence of the O enhancement in 
the solid phase. Consequently, the chemical abundances 
(an their distribution throughout the fluid layer) of both 
C and O are expected to be significantly different from 
those of the previous model. The Brunt- Vaisala frequency 
corresponding to a CO white dwarf model in which the 
chemical profile has been computed self-consistenly with 
the predictions of chemical rchomogenization (Fig. |2Jl is 
shown in Fig. Note that at this stage, the CO white 
dwarf models have a significant percentage of its interior 
crystallized (55%), in such a way that the chemical step 
produced by core overshooting during the pre-white dwarf 
evolution has been wiped out by the ongoing chemical re- 
distribution. Because the external chemical transition re- 
gions have been smoothed by diffusion, the shape of the 
Brunt- Vaisala frequency has not any abrupt feature ca- 
pable to strongly pertubate the period spectrum of the 
model. Hence, the period spacing distribution should be 
more or less uniform, at least in the asymptotic limit of 
large overtones (long periods). 

Now, we examine the Brunt-Vaisala frequency charac- 
terizing a typical ONe white dwarf model at the ZZ Ceti 
stage. In Fig. |7| we plot the profile of N"^ for the same 
model analyzed in Fig. 01 Clearly, the dominant feature 
is that located at —2 <^ log(l — M^/Mf,) <^ —1 (see the 
Ledoux term B in the inset of the figure). As we shall see 
in the next section, this feature is responsible for strong 
non-uniformities in the period spectrum of this model. 

3.2. Pulsational spectrum 

We have computed adiabatic, nonradial, spheroidal 
(7(gravity)-modes with £ = 2 covering the period range 
of pulsations observed in ZZ Ceti stars. We have not con- 
sidered torsional modes, since these modes, characterized 
by very short periods — of up to 20 s (Montgomery & 
Winget 1999) — have not been detected in ZZ Ceti stars 
so far. 

We begin by examining Fig. |H1 where the asymptotic 
period spacing as a function of the effective temperature 
is shown for the two CO and the ONe white dwarf mod- 
els previously discussed, and for two additional models 
in which crystallization has been completely disregarded. 
These two additional models are only computed for the 
sake of comparison and in order to exemplify the effects 
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Fig. 7. Same as Fig.lHl but for the ONe white dwarf model 
analyzed in Fig. ^ Dashed Hnes correspond to the case 
in which the Ledoux term B is neglected in the region 
-2 ^ log(l - M./NU) <^ -1 (see Sect. 3.2 for additional 
details). 

of crystaUization in spite of the fact that they are not 
physically sound. The thick lines correspond to our fidu- 
cial models. That is, the thick solid line corresponds to an 
ONe model in which crystallization has been taken into 
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Fig. 8. Asymptotic period spacing, AP(^^2)j as a func- 
tion of the effective temperature, for the CO and ONe 
evolutionary cooling sequences (thick and thin curves, re- 
spectively). See text for details. 



account, whereas the thick dashed line corresponds to the 
CO white dwarf model in which both crystallization and 
phase separation were taken into account. The thin dashed 
line corresponds to a CO model in which crystallization 
but no phase separation was considered. On its hand, the 
thin dotted line corresponds to a CO model in which nei- 
ther crystallization nor phase separation were accounted 
for. Finally, the thin solid hne indicates the results for an 
ONe model in which crystallization was completely dis- 
regarded. The results presented in Fig. |H1 deserve several 
comments. Firstly, for all the cases considered AP^=2 is 
an increasing function of Tgfi, as expected from the fact 
that the Brunt- Vaisala frequency decreases as the white 
dwarf cools down — Xt decreases as degeneracy increases, 
see Eq. Secondly, it is noteworthy that in the case in 
which crystallization is neglected APg=2 is smaller for the 
ONe models than for the CO ones. This trend is under- 
stood in terms of N"^ having notably larger values at the 
core region in the ONe models compared to the CO mod- 
els. Thus, the integral of Eq. turns out to be larger in 
the case of ONe models and, consequently, the asymptotic 
period spacing is smaller. 

In contrast, for the cases in which crystallization has 
been considered, the values of APg=2 of the ONe models 
are always larger than those of the CO ones. In part, this 
can be understood on the fact that the ONe model is char- 
acterized by a larger crystallized mass fraction than the 
CO one for a fixed Toff- This in turn implies that the inte- 
gral of N/ r in Eq. is smaller and the resulting asymp- 
totic period spacing in the ONe model is larger. With re- 
gard to the CO models, we note that, when chemical reho- 
mogenization is considered the asymptotic period spacing 
is slightly greater when compared with the situation in 
which rehomogenization is neglected. This can be under- 
stood on the basis that, in the case in which rehomogeniza- 
tion is not allowed to operate, the integral in Eq. has 
an extra contribution from the overshoot-induced step via 
the strong peak at log(l — Mi/A/») « —0.7 in the Brunt- 
Vaisala frequency (see Fig. 

We now examine the forward period spacing, defined 
as APk = Pk+i — Pk (being k the radial overtone). This 
quantity, which is particularly sensitive to the details of 
internal chemical structure of the white dwarf models, is 
usually employed to infer the mode trapping properties of 
ZZ Ceti stars. In addition, the theoretical period spacing 
can be directly compared with the observed one in a given 
star, whenever a sufficient number of consecutive periods 
(with the same harmonic degree €} can be measured. Here 
we restrict ourselves to present pulsational results corre- 
sponding to the same models analyzed in Figs. [S] and 
H We recall that these models are characterized by an ef- 
fective temperature of « 11800 K, which is representative 
of ZZ Ceti stars. Fig. El shows the i = 2 period spacing 
distribution (top panels) and the kinetic energy of oscilla- 
tion (bottom panels) as a function of the period for a CO 
white dwarf in which crystallization has been ignored (left 
panels), for a CO white dwarf in which crystallization was 
taken into account but no phase separation was consid- 
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Fig. 9. The forward period spacing (upper panels) and the kinetic energy (lower panels) in terms of the periods of 
£ = 2 pulsation modes. The left panels correspond to a CO white dwarf model in which crystallization has not been 
taken into account. The central panels show the same quantities for a model at the same Teff in which crystallization 
but no phase separation have been considered. Finally, the right panels show the case in which both crystallization 
and phase separation have been accounted for. 



ered (middle panels) and for the case in which both crys- 
tallization and phase separation were properly accounted 
for (right panels). 

Clearly, when crystallization is considered the period 
spacing distribution is much less featured. A similar be- 
havior has been reported by Montgomery & Winget (1999) 
in the frame of calculations in which the crystallized mass 
fraction is considered as a free parameter. When crystal- 
lization is ignored, the internal boundary condition re- 
mains fixed at the center of the star. The resulting eigen- 
spectrum is composed by some core-confined modes, some 
modes trapped in the outer layers, and the remainder nor- 
mal modes (Corsico et al. 2002; Althaus et al. 2003). In 
this case the modes nearly preserve their character during 
the evolution. In contrast, when crystallization is consid- 
ered (middle panels) the inner boundary condition moves 
together with the crystallization front, and the eigenfunc- 
tions of all modes are pulled out to the surface of the 
model. Under these circumstances, the resonance condi- 
tions determining which modes are trapped or confined 
are strongly modified as the degree of crystallization in- 
creases. This is clearly shown in the top middle panel of 
Fig. ini which corresponds to a model in which the crys- 



tallization front is located at Mc — 0.63 Af*. Note that 
for the high- and intermediate-overtone modes the strong 
minima (associated with core-trapped modes) character- 
izing the APfe distribution of the model for which crys- 
tallization was ignored (left panel), are not present in the 
case in which crystallization is considered (middle panel). 
A look to the lower panels of Fig. |51 shows as well that 
when crystallization is taken into account the differences 
between the kinetic energies of trapped, normal and con- 
fined modes are smaller. Another interesting point is that 
the mean period spacing increases when crystallization 
is considered, in agreement with the predictions of the 
asymptotic theory of pulsations and with the results of 
Montgomery & Winget (1999). 

We now examine the more realistic case in which phase 
separation is included in addition to crystallization (right 
panels). In view of the above considerations, we expect a 
smooth period spacing distribution. We note that the pe- 
riod spacing distribution (upper right panel) has a fairly 
clean structure for periods longer than « 300 s. In fact, 
the forward period spacing for modes with fc > 12 becomes 
very close to 20.6 s, the value predicted by the asymptotic 
theory, given by Eq. Q). As already discussed, this is due 
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mostly to the absence of any abrupt feature in the profile 
of the Brunt- Vaisala frequency. In fact, the model is char- 
acterized by a crystallized mass fraction of « 0.55, which 
is large enough for the ongoing chemical rehomogeniza- 
tion processes to erase the overshoot-induced step in the 
chemical profile (see Figs. Q] and [2Jl. 

As previously stated, the major aim of this work is to 
explore the possibility of using white dwarf asteroseismol- 
ogy to distinguish massive CO white dwarfs from those 
having ONe cores. In the right panels of Fig. ^| we show 
the distributions of AP^ and ii^kin corresponding to our 
fiducial ONe model at Teff ~ 11800 K (which corresponds 
to a crystallized mass Mc f» 0.90 M*). For the sake of com- 
pleteness in the left panels of this figure we also show the 
same distributions for the case in which crystallization 
was ignored. We note that, irrespective of the assump- 
tions regarding crystallization for both types of models 
(CO or ONe), the period spacing diagrams of ONe mod- 
els are notably different from those of the CO models. In 
particular, for the case in which crystallization was dis- 
regarded (left panels), the period spacing diagram of the 
ONe model shows strong and abundant features linked 
with mode trapping and confining, mostly driven by the 
structure at -2 ^ log(l - M^/M*) ^ -1 in the Brunt- 
Vaisala frequency profile (see Fig. CJ). When crystalliza- 
tion is taken into account, the period spacing distribution 
of the ONe model still has a complex structure but with 
a notably lower strength (right upper panel of Fig. ITUj) . 
In particular, we note the presence of three minima at 
« 300, 600 and 900 sec. However, the comparison of the 
period spacing diagrams in the right panels of Figs. El and 
EH indicates that two models with the same stellar mass 
and effective temperature but differing in the core compo- 
sition (one having a core rich in carbon and oxygen, and 
the other mainly composed by oxygen and neon) should 
be characterized by a quite different pulsational spectrum. 

An additional difference in the pulsational properties 
between both types of models is found when the kinetic 
energies of the modes are compared. In fact, all the modes 
of the ONe model have lower energies than those of their 
counterparts of the CO model. Two facts help in clarify- 
ing this issue. We recall that Skin is proportional to the 
density and to the squared radial and horizontal eigen- 
functions. Firstly, the propagation region in which the 
eigenmodes are allowed to oscillate have different sizes. 
Specifically, for the ONe white dwarf at Toff ~ 11800 K, 
we find that 90% of the stellar mass has already crystal- 
lized. Hence, in this high-density region nonradial g-modes 
are excluded. Therefore, the modes are capable to propa- 
gate only in the remaining external, low-density region. In 
the case of the CO model at the same Toff, the crystallized 
mass fraction is aproximately 0.55, and so the modes can 
propagate in regions characterized by considerably high 
densities. Secondly, we have found that for the case of 
ONe models the amplitudes of the eigenfunctions are con- 
siderably smaller than those of their counterparts in the 
CO models. As a final remark, we note that the mean pe- 
riod spacing of the pulsation modes corresponding to ONe 
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Fig. 11. Same as the upper-right panel of Fig. ^| but 
for the case in which the effect of the strong step in the 
chemical profile at —2 ^ log(l — M-^/M^,) <^ —1 (see Fig. 
0J) on the shape of the Brunt- Vaisala frequency has been 
neglected by setting _B = in that region (see Fig. d). 



white dwarfs is slightly larger than that of the CO mod- 
els. This is again in agreement with the predictions of the 
asymptotic theory, as shown in Fig. |S1 

In closing, we stress that the main differences in the 
pulsational spectrum between ONe and CO models are 
due to the presence of a strong step in the core chemical 
profile of the ONe model. As mentioned, this step marks 
the outer edge of the mixing zone induced by the negative 
molecular weight in the deeper layers of the ONe white 
dwarf progenitor. However, it could be possible that ad- 
ditional mixing episodes beyond this external border may 
wipe out the chemical step. This being the case, we would 
expect the period pattern to be modified to a consider- 
able extent. This is indeed borne out by Fig. 1111 in which 
the period spacing is computed under the assumption that 
B = in the region at -2 ^ log(l - M^/M^) ^ -1 and 
the resulting run of the Brunt- Vaisala frequency is that 
displayed with dashed line in Fig. [71 We are aware that 
this procedure docs not eliminate completely the effects 
of the strong chemical feature on the period spectrum, 
because the change in density associated with this chem- 
ical interface is still present in the computation of N (see 
Eq. 2). However, we consider this approximation reliable 
enough for our purposes here. Note that the forward pe- 
riod spacing becomes much less featured as compared with 
the case in which the Brunt- Vaisala frequency is computed 
self-consistently (right panel of Fig. I10|) . and instead, it 
is reminiscent of the results obtained for the CO model 
(see right panel of Fig. O. This aspect would render the 
current asteroseismological techniques much less capable 
to infer the actual core composition of the massive white 
dwarfs. We note, however, that the average period spac- 
ing remains virtually unchanged even in presence of such 
extra mixing episodes. 
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Fig. 10. Same as Fig. but for the case of an ONe white dwarf. 
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4. Conclusions 

In this work we have computed the evolution of mas- 
sive white dwarfs with hydrogen-rich envelopes and 
carbon/oxygen and oxygen/neon cores with the major 
aim of comparing their adiabatic pulsational properties. 
Particular attention was given to the formation of the core 
chemical profile during the pre-white dwarf evolutionary 
stages and to the evolution of the chemical abundance 
distribution caused by element diffusion during the white 
dwarf regime. In addition, the chemical rehomogenization 
induced by phase separation upon crystallization as well as 
the effects of crystallization on the oscillation eigenmodes 
have been fully taken into account. 

Our results show that because the chemical inter- 
faces in the envelope are markedly smoothed by diffu- 
sion, the chemical structure of the core mostly deter- 
mines the structure of the period spectrum. This result 
reinforces the conclusions arrived at in Corsico et al. 
(2001), (2002) and Althaus et al. (2003) about the role 
of time-dependent element diffusion in white dwarf pul- 
sation calculations. Moreover, we find that, in contrast 
with their carbon/oxygen core counterparts, oxygcn/ncon 
white dwarfs are characterized by strong non-uniformities 
in their period spacing distribution. The lack of a featured 
period spacing distribution characterizing white dwarfs 
with carbon/oxygen cores arises from the fact that the 



mixing episode induced by phase separation as crystal- 
lization proceeds completely wipes out any feature in the 
shape of the innermost chemical profile built up during the 
evolutionary stages prior to the white dwarf formation. In 
addition, we found that the average period spacing char- 
acterizing our oxygen/neon models is appreciably larger 
than that of the carbon/oxygen counterparts having the 
same stellar mass and effective temperature. This is true 
even in the presence of extra mixing episodes that could 
eventually remove the sharp abundance discontinuity left 
in the oxyen n/neon models. Finally, we also find a quite 
different kinetic energy spectrum for both types of models. 

Hence, it turns out that the asteroseismological tech- 
niques usually employed in variable white dwarf studies 
are, in principle, a powerful tool for distinguishing mas- 
sive carbon/oxygen white dwarfs from those having oxy- 
gen/neon cores. In this regard, the most massive ZZ Ceti 
star presently known, BPM 37093, is of particular interest 
since the value of its stellar mass of 1.05 Mq places it 
close to the theoretical lower limit for a oxygen/neon white 
dwarf to be formed. Nonetheless, the scarcity of detected 
periods (eight periods) in the light curve of BPM 37093 
renders conclusive asterosesimological inferences in this 
case very difficult. However, since the average period spac- 
ing is larger for the oxygen/neon models than for the car- 
bon/oxygen ones, then a comparison with the observed 
value would help to get a clearer insight into the actual 
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core composition of this massive white dwarf. We feel that 
this aspect would deserve to be explored in the frame of 
a set of white dwarf models with various stellar masses. 
Ongoing attempts to constrain the fundamental proper- 
ties of BPM 37093 using asteroseismological data would 
allow us to shed new light on the scenarios leading to the 
formation of single massive white dwarfs and on the na- 
ture of their progenitors. 
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